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We reconsider the recent observation by the D0 experiment of a sizable like-sign dimuon charge 
asymmetry, highlighting that it could be affected by CP- violating New Physics contributions not only 
in B<i- and l? s -meson mixings, but also in semileptonic decays of b and c quarks producing muons. 
The D0 measurement could be reconciled with the Standard Model expectations for neutral-meson 
mixings, provided that the CP asymmetry in semileptonic b (c) decays reaches 0.3 % (1%). Such 
effects, which lie within the available (rather loose) experimental bounds, would be clear indications 
of New Physics and should be investigated experimentally. 



1 Introduction 

The recent observation by the D0 experiment of a sizable like-sign dimuon charge asymmetry [1-3] 
has attracted tremendous theoretical interest. If interpreted as originating from mixing- induced CP 
violation in semileptonic b quark decays, the observed value of Ref. [3] 

A b sl = (-0.787 ±0.172 ±0.093)%, (1) 

is significantly larger than expected within the standard model (SM) .A^ = (— 3.96^q'q4) • 10 -4 [4] 
(see also Refs. [5]). On the other hand, recent measurements of CP violation in wrong sign semilep- 
tonic B s decays from the LHCb and D0 collaborations [6, 7] 

a* A = (-1.08 ± 0.72 ± 0.17)% [DO] , a s A = (-0.24 ± 0.54 ± 0.33)% [LHCb] , (2) 

are in agreement with SM expectations. Furthermore, recent precise determinations of the CP asym- 
metry in the B s —> J/ip(j> decay [8], which are consistent with SM predictions, already severely 
constrain a possible interpretation of the value in Eq. (1) in terms of non-standard CP-violating con- 
tributions to B s mixing [4]. While the large observed could still originate from new CP-violating 
effects in B^ mixing, this would nonetheless require sizable correlated non-standard contributions 
also to the absorptive mixing amplitude [4] (see also Ref. [9]). 

In view of this intriguing situation, it is important to revisit theoretical assumptions underlying 
the interpretation of the D0 measurement. Since the D0 result indicates a significant discrepancy 
with the SM concerning a tiny CP asymmetry, one should reconsider other possible New Physics 
(NP) sources which could contribute the observed CP violation, but are generally assumed to be 
negligible within the SM. Indeed one can imagine other sources of CP violation contributing to 
the dimuon charge asymmetry. Given the sizable charm contributions in both the inclusive muon 
and same-sign dimuon data samples, an attractive possibility is represented by CP violation in D° 
mixing. However, the D° oscillation probability has been measured and its smallness renders any 
such contribution negligible. 

On the other hand the D0 analysis also assumes that the semileptonic b and c quark decays 
conserve CP. In this letter we explore the consequences of relaxing this assumption. Direct CP 



violation in semileptonic decays of b or c quarks is presently only poorly constrained experimentally. 
The related inclusive semileptonic CP asymmetries are expected to be tiny in the SM and so offer 
interesting venues for contributions of NP. This motivates us to investigate whether the observed 
dimuon asymmetry in Eq. (1) could originate from CP violation in semileptonic b and c quark 
decays. 

The rest of the paper is structured as follows: in Sec. 2 we recall the basic elements of the D0 
analysis needed for the extraction of the CP asymmetries from their inclusive muon and same-sign 
dimuon data samples. Sec. 3 is devoted to the derivation of the relevant observables in presence of 
CP violation in semileptonic b and c quark decays and the reinterpretation of the D0 measurement 
in terms of the related CP asymmetries. In Sec. 4 we discuss other existing bounds on CP violation 
in these decays, and we give their expectations within the SM in Sec. 5. Finally, we conclude in 
Sec. 6. 



2 Elements of the experimental analysis 

We recall the basic ingredients in the determination of from the measurement of the dimuon 
charge asymmetry according to Ref. [2] and applicable to the most recent update [3] (we will use the 
values quoted in the latter reference for our numerical analysis). Experimentally, both the like-sign 
dimuon charge asymmetry A and the inclusive muon charge asymmetry a are measured by counting 

N++ - N~ _n+ -rT 

A ~ N++ + N— ' a ~ n + + n- ' (3) 

where N ++ ~~ (n + '~) denote the number of events with two muons (one muon) with a given charge 
passing the kinematic selections. 

There are various detector and material-related processes contributing to these asymmetries. The 
reconstructed muons are classified into two categories: "short" (or S) including muons from weak 
decays of b,c,r and from electromagnetic decays of short-lived mesons (0, uj, rj, p°), and "long" (or 
L) coming from decays of charged kaons and pions as well as from particle misidentification. One 
can separate the contribution from short muons as to the inclusive muon charge asymmetry 

a = fs(as + o~) + fKCLK + /ttOtt + fpdp , (4) 

where 5 is the charged asymmetry related to muon detection and identification, and for the back- 
ground L processes, fx (f-K, f P ) the fraction of muons from charged kaon decays (charged pion 
decays, misidentifications) , and clk (a n ,a p ) the corresponding charge asymmetry. These quantities 
are directly measured from experiment. 

One has also a similar expression for the dimuon charge asymmetry 

A = F SS A S + (F bkg - 2F LL )a s + (2 - F bkg )A + F K A K + F„A W + F p A p , (5) 

where A is the detection asymmetry, Fx t7T ,p and Ax^^ are the fractions and asymmetries related to 
the various background L processes. 

Finally, D0 also consider the combination 

A' = A - aa = F SS A S + (F bkg - 2F LL - af s )a s + ..., (6) 

where a = 0.959 [2] or a = 0.89 [3] is tuned to reduce the background uncertainties (denoted by the 
ellipses) , leading to a precise extraction of a linear combination of as and A s . 

Various processes (T^ ) producing fj,^ can contribute to the asymmetries as and ^ with relative 
weights (wi) as determined from Monte Carlo simulations, see Table 1. We start with the probabilities 
^b{c) f° r an initial 6(c) quark in bb production (prompt cc production, without b hadrons) to produce a 
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Table 1: Processes contributing to the inclusive muon and like-sign dimuon samples, adapted from 
Refs. [2,3] to include CP violation in Bd,B s mixings (A^) and in semileptonic decays (A dix , A c dh ). 
Numerical values of the weights (normalized to w\ = 1) used in the analysis are W2 = 0.096 ± 0.012, 
7ij 3 = 0.064 ±0.006, w 4 = w 4 a + w 4b + w 4 c = 0.021 ±0.001, w 5 = 0.013 ±0.002, u> 6 = 0.675 ±0.101 [3], 
and the mean B ds mixing probability Xo = 0.147 ± 0.011 [2] or Xo = 0.1259 ± 0.0042 [3]. 
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H^, respectively. In the same processes, P b /A are t nen the corresponding probabilities to also produce 
a from the accompanying 6(c) quark. Finally we have the probabilities P^ M = P^m to produce 
ji^ from short-lived meson (or SLM: 77, ui, p°, <fr, J/ip, if)') in events with no b or c hadrons (we assume 
these processes to conserve CP). We should mention that this table is adapted from Refs. [2,3] by 
distinguishing the CP-conjugate processes in order to include the CP- violating effects discussed in 
the present letter. 

Assuming all the processes producing muons to be independent 1 , we can now construct both 
asymmetries as: 

= E 9=fe , c , S LM[(^ + + p Q + ) - cgr + gr)j (7) 

aS T, q=b , c ,SLMl(P? + Pf) + (P? + ^)V 
and _^ 

52q=b,c,SLM[(Pg~ P g) ~ (Pq ' P q 

Y^ q =b,c,SLM[(P<t " Pq) + (Pq ' Pq 

Let us finally mention that another observable is introduced in Ref. [3] by considering A' s restricted 
to a sample of dimuons with a large enough muon impact parameter. We refrain from studying 
in detail this observable because of our lack of knowledge concerning the experimental inputs and 
correlations required, but we highlight that it could be analyzed along the same lines as what we 
present here for A' s . 

3 The role of the semileptonic CP asymmetries 



. ^g=6,c,!>LMLV'9 q ' \ q q Ji ,„s 
A S = ^ „ n + ^ i 77^ fS-M • ( 8 ) 



In the absence of CP violation in semileptonic 6 and c quark decays, 05 and As can be directly 

ib 

si" 



expressed in terms of A**,, a linear combination of the wrong-sign semileptonic flavor specific asym- 



metries of the P>d s mesons measuring CP- violation in their respective mixings 



2. 



A sl - f d a sl + f s a sl , a sl - ^ - — } - ^ , (9) 

with fd and f s the fractions of B& and B s mesons contributing to the asymmetry, which depend on 
the experimental setting. The values used by the D0 collaboration are either taken from averages 
at Tevatron [2], or at LEP [3]. 

In the presence of CP violation in inclusive semileptonic 6 or c decays, one must define two 
additional asymmetries A dh and A div : 

b T(b^fi-x)-r(b^fi+x) r(c->/z-jr)-r( c 

^ F(6^ fI -X) + T(b^ n+X)' dir T(c^ fi-X) + T{c^ fi+X) ' { ' 

Then the most general expressions for the probabilities P^ , P^ read 

P+ oc w lb (l + Al l -A b dh ) + w 2 a(l-A c dir ) + (w 3 + w 5 )(l-A c dir )/2 + w± a /2, (11a) 

P- OC ^i a (l+^+W26(l + ^sl + ^dir) + (^3 + ^5)(l + ^dir)/ 2 + W 2 > ( llb ) 

P+ oc w la (l-A dir ) + W2 b (l-A^-A c diT ) + (w 3 + w 5 )(l-A c diT )/2 + w^/2, (11c) 

P b ~ oc w u (l-Al l + A dir ) + w 2a (l + A c dir ) + (w 3 + w 5 )(l + A c dir )/2 + w 4a /2, (lid) 

P+ oc w 6 (l-A c div ) + w 4b /2, (lie) 

P- oc w 6 (l + A c dir ) + w 4b /2. (llf) 



1 The correctness of this approximation was verified numerically using Monte Carlo generated data samples in Ref. [2] 
for the case of CP violation in Bd, s mixing. 

2 As mentioned in the Introduction, possible contributions due to CP violation in D° mixing are extremely suppressed 
by the small mixing probability of D° mesons as measured experimentally. 
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CP asymmetry (%) 


a 


A 


A' 


Ki (^dir = 0,^ ir = 0) [3] 


-1.04 ±1.30 ±2.31 


-0.808 ±0.202 ±0.222 


-0.787 ±0.172 ±0.093 


< (4 6 l = 4 6 l SM ^Sir = 0) 


0.11 ±0.15 ±0.28 


0.26 ±0.07 ±0.07 


0.28 ±0.06 ±0.04 


^dir (4 6 1 = ^1 SM ^ir = 0) 


0.13 ±0.17 ±0.31 


0.69 ±0.18 ±0.21 


0.93 ±0.21 ±0.20 



Table 2: CP asymmetries (in %) needed in meson mixing (first row) or semileptonic decays (second 
and third rows) in order to obtain the measured values of a, A, A' from Ref. [3]. 



Furthermore, since semileptonic charm decay contributions to wrong-sign muons P c and P^ are 
suppressed by the small D° mixing probability, we have simply 

P~ =P+ ocuW2, (12) 

whereas the short-lived meson decays provide 

^SLM = P SLM « w 4c/2 . (13) 

We can then use Eqs. (7) and (8) to express the three observables related to S muon production 
a s (A b ah A c dk ,A b dir ), A s (A b sl ,A c dir ,A b dir ), A' s (A b sl , A c dir , A b dir ) = F ss A s + (F bkg -2F LL -af s )as , (14) 

in terms of the three asymmetries of interest. 

Let us emphasize that in principle, one needs to distinguish the weights corresponding to the 
three different contributions from short-lived mesons T^ a ^b^ c . Unfortunately, D0 does not provide 
the three probabilities separately. By varying them in the intervals W4 C £ [0,104], w^b £ [0, 104 — W4 C ] 
and w^a £ [0,104 — ^4c — Wib], we have checked however that the associated additional uncertainty 
in relating A b v A dh and A dh to as and As is negligible (for defmiteness, we present the results 
corresponding to = w^ c = in the following). In the limiting case where A dir = A dir = and 
w 4b = u>4 C = 0, the resulting expressions coincide with the corresponding expressions in Ref. [3]. 

Unfortunately, lacking complete information on the correlations among the various inputs entering 
the D0 measurement, we cannot extract the values of A dil and/or A dir from the measurements of 
a, A (and A') in a trustable way. However, using as, As and A' s , we are in the position to at least 
estimate the values of A dh or A b dir assuming A b \ as predicted in the SM 3 . We do this by looking for 
solutions of the following equations 

asiA^, A c dir , 0) = a 5 « a , 0, 0) , (15a) 
A 5 (^ SM ,^ ir ,0) = A s (A b slA ,0,0) , (15b) 
A' 5 (^ SM ,^ ir ,0) = A s (A b l Af , 0,0) , (15c) 

and similarly for the case of nonzero A dh instead of ^ir- Here A^ a ,A b l A , A b { A , are values of A^ as 
extracted from a, A or A' in Ref. [3]. The results are collected in Table 2. The uncertainties quoted are 
dominated by those coming from A b ^ a , A b \ A , A b y A ,. These implicitly include the uncertainties (and 
correlations) from the other parameters entering the D0 analysis. We furthermore include explicit 
contributions from additional input parameters entering equations (15a)-(15c) to the systematical 
uncertainties (without the proper knowledge of their correlations with A b sl a , A b sX A , A b sl A , , this can 
only be considered a very rough estimate of the potential size of their effects). We note that these 

3 We consider the SM result from Ref. [4], which is larger in size than the SM result quoted in Refs. [2,3], and thus 
our results can be considered more conservative concerning the required size of CP-violation in semileptonic b or c 
decays. 
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additional contributions only significantly affect the extraction of A% from A', where they almost 
double the systematic error budget. On the other hand, the uncertainties coming from the SM 
prediction for A^ are completely subleading and thus not quoted. Finally, we have checked that 
consistent results (but with slightly larger errors) are obtained when using inputs and asymmetry 
measurements from the previous D0 analysis [1,2]. 

One notices that CP asymmetries in inclusive semileptonic b (c) decays below the 0.3%(1%) level 
are required to explain the D0 measurement. We also note that the values of A dh as extracted 
from a and A are well consistent at OAa level (assuming Gaussian uncertainties and in absence of 
correlations), while there is a slight 1.2a difference between the values of A^ ir extracted the same way. 
Although the values extracted from A' are even bigger, they are expected to be highly correlated 
with the ones from the other two observables, and we do not attempt to assign a statistical meaning 
to these differences. 

One can imagine that both A b Air and A^ ir may differ from zero. It is then useful to compute the 
dependence of as, As and A' s on the three CP asymmetries (to first order), yielding: 

a s = 4^(0.061 ± 0.004) +^ ir (-0.535± 0.028) + A| ir (-0.454 ± 0.028), (16a) 
A s = A^(0.474 ± 0.023) +^ ir (-1.421± 0.024) + ^(-0.527 ± 0.025), (16b) 
A's = 4^(0.312 ± 0.023) +^ ir (-0.849± 0.061) + A| ir (-0.250 ± 0.038), (16c) 

where the relevant systematical and statistical uncertainties have been combined in quadrature. We 
see that the observables as, As, A' s exhibit similar sensitivities to the three types of CP violation. 
This clearly indicates that the interpretation of these quantities in terms of neutral-meson mixing 
requires a further check of the absence of CP violation in decays at a similar level to the uncertainties 
quoted for A^. 

Assuming the SM value of A^,, the experimental values of as and As set constraints in the 
(A^ ir , A^ ir ) plane, as illustrated in Fig. 1. One can see that the sensitivity of the observables to A h dlr 
is larger than that to Aji-, explaining that a larger asymmetry in charm is required to reproduce the 
D0 value for the dimuon asymmetry. Since our analysis does not include all the relevant correlations, 
we do not attempt at combining the two constraints statistically, even though this could be done 
easily by the D0 collaboration. 

4 Other Constraints on Direct CP violation in semileptonic b and 

c decays 

At present and to the best of our knowledge, there exist no direct experimental bound on the 
CP-violation in semileptonic decays of either b or c quarks [10, 11]. We can nonetheless form a 
poor man's estimate of the current constraints by considering the uncertainties on the CP averaged 
(semi)leptonic branching fractions (estimating in this way the amount of CP-violation in semileptonic 
decays that could be hidden by the systematics of the measurement). In particular, among the various 
measured semileptonic D + , D° and D s decays, the CP averaged branching fraction of D + — > K^p^v^ 
(which accounts for more than half the inclusive semileptonic D + branching fraction) is currently 
most precisely known with 6.5% relative uncertainty [12]. Based on this and considering also the 
experimental bound on the CP asymmetries in leptonic — > y^v decays [13] we observe that at 
the current level of uncertainty, CP-violation in semileptonic D decays at or below 6% is still viable. 

Concerning leptonic and semileptonic decays of B mesons, we find that the current level of 
uncertainty of the CP averaged inclusive B -4 X c lv branching fractions [10] in principle still allows 
for CP violation in semileptonic B decays at the level of around 3%. However in this case, additional 
constraints can be derived using the measured flavor specific semileptonic charge asymmetries defined 
in Eq. (9) [11], which yields for the B& meson 

a d sl = (-0.05 ± 0.56)% , (17) 
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Figure 1: The A c Air and A b dir semileptonic decay asymmetries needed in order to reproduce the 
measured values of inclusive semileptonic and same-sign dimuon asymmetries a (in thick contours) 
and A (in thin contours) (results for A 1 closely resemble those for A) using the SM predicted value 
for the inclusive wrong-sign semileptonic asymmetry A b , (the weights and parameters of Ref. [3] have 
been used). The la (2a) bands are marked with full (dashed) contours. 



whereas the situation for a s sl is recalled in Eq. (2). These measurements are generally based on 
specific charmed decay modes associated with muons, so that they can be affected by CP violation 
in semileptonic b decays, but not in c decays, contrary to the D0 measurement. 

In reinterpreting these results in terms of CP-violating asymmetries in inclusive semileptonic b 
decays at D0 , we first identify |o^| > lAiirl an ^ \a s s i\ ^ l-^dirl' wriere the semileptonic asymmetries 
^■dir are defined as in Eq. (10), but now refer to the relevant decaying B q mesons. Next we need to 
sum over the relative production fractions of the various b hadrons contributing to the semileptonic 
event samples: 

A b dil = f(B u )A B d - + f(B d )A^ + f(B s )A^ + . . . , (18) 

where the ellipses denote neglected smaller b hadronic state contributions. As a first approximation 
for f(B q ), we neglect different lifetime effects and use the measured unbiased 6-hadron fractions as 
measured by various experiments at high energies /(Bj) = f(B u ) = 0.401 ± 0.007 and f(B s ) = 
0.107 ± 0.005 [11]. Neglecting possible correlations between the various inputs and combining all 
uncertainties in quadrature we obtain a bound of |^4 dir | < 1.2% which is safely above what is required 
for explaining the D0 result. 



5 SM Expectations for A b dil . and A c Ah: 

Finally let us briefly comment on the expect size of A, ir and A dil within the SM. Direct CP violation 
in decays requires the presence of (at least) two interfering decay amplitudes (we will denote them 
as At = \At\ expi((pT + <5t) and Al = \Al\ expi((pL + 5l)) with different weak (<^t,l) and strong 
($t,l) phases. Denoting 1Z = \Al/At\, = 4>l~ 4>t and A<5 = 5l — 5t, the related CP asymmetry 
can then be written as 

2ft sin sin A(/> 
dir ~ 1 + 211 cos A5cosA0 + ft 2 ' ( ' 



7 



The dominant tree-level SM contributions to semileptonic transitions (below the W scale) are 
described by the relevant effective weak Hamiltonian 



si _ 
eft — 



AG F 



E E W£Ml-75)£][W-75N+h.c.. 



(20) 



y/2 



U =u,c D=d,s,b 



The required additional amplitudes can be generated at one-loop via time ordered correlators of 
H^ff with the effective weak Hamiltonian describing non-leptonic b — > cud and b — > cus decays 
J d 4 xT{%"g(0), %|g(x)}. The presence of light quarks in the loop provides a source for the strong 
phase difference, while the weak phase difference is encoded in the relevant CKM matrix elements. 
However, being higher order in Gf ~ ^/ v ew these effects are expected to be severely suppressed 
leading to 1Z <C 1. Expanding Eq. (19) to linear order in 7Z, assuming sin = 0(1) and using naive 
dimensional analysis, we can estimate the resulting SM contributions to A^ iv and A^ iv as 



We observe that the sizes of A^ ir or A^ iv required to accommodate the D0 result are orders of 
magnitude larger than the above SM expectations for these quantities and their confirmation would 
thus constitute a clear indication of NP. 

6 Conclusions 

We have reconsidered the recent measurement made by the D0 collaboration of a like-sign dimuon 
asymmetry [1-3]. Since this measurement disagrees significantly with the SM prediction for CP 
violation in Bd and i^-meson mixing, we have reassessed some of the underlying assumptions of the 
analysis, allowing for NP effects violating CP not only in mixing, but also in decays. 

We have shown that the D0 result can be made compatible with the SM expectations for CP 
violation in Bd s mixing, provided that non-SM contributions are introduced either in 

• CP violation in b semileptonic decays (-A^ir ~ 0.3%), which is currently allowed experimentally, 
though close to the uncertainties quoted for the individual semileptonic Bd and B s asymmetries. 

• CP violation in c semileptonic decays (Aj iT ~ 1%), which is currently allowed experimentally. 

As indicated in Fig. 1, one could also consider the presence of both effects. We discussed briefly 
the size of these effects in the SM, and showed that they are much smaller than what is required to 
explain the D0 result. 

At this stage, the possibility of a CP violation in semileptonic charm decays (at the percent level) 
is actually quite intriguing and one could speculate about the connection of such effect with the 
significant CP violation observed in hadronic charmed decays [14,15]. In any case, a CP-violating 
contribution to charm semileptonic decays would allow the D0 measurement of the like-sign dimuon 
asymmetry to differ from the SM value, but it would not affect the measurements of or, based on 
specific decay channels like B s — > D s fiX. 

Our analysis is obviously very naive as far as experimental uncertainties and correlations are 
concerned. The numbers provided in this short note are only indicative, but we firmly hope that they 
will incite experimentalists to revisit the D0 analysis and related studies to include and constrain CP 
violation in semileptonic b and c decays. Such cross-checks would be particularly useful to improve 
our understanding of neutral- meson mixing and its potential for NP searches. 
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